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Abstract 

We have developed a theory of three-pulse coherent control of photochemical processes. It is 
based on adiabatic passage and quantum coherence and interference attributed to the lower-lying 
dissociation continuum and excited upper discrete states, which are otherwise not connected to the 
ground state by one-photon transitions. Novel opportunities offered by the proposed scheme are 
demonstrated through extensive numerical simulations with the aid of a model relevant to typical 
experiments. The opportunities for manipulating the distribution of the population among discrete 
and continuous states with any necessary ratio by the end of the pulses are demonstrated. 
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I. INTRODUCTION 



Laser control of chemical reactions and other dissipation processes like photodissociation 
of molecules and photoionization of atoms has recently gained increa sing interest in the 
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laser and chemical communities 

lisl . Q]. In many cases, such control is based on the interference of different quantum 
pathways. While the main stream of theoretical and experimental work on coherent quantum 
control concerns laser-induced transitions between discrete states, coherence and interference 
effects attributed to dissociation and ionization energy continua are of great importance, 
because of various applications in atomic and molecular physics and photochemistry. The 
properties of the continuum of quantum states such as observed in the ionization of an atom 
or dissociation of a molecule have attracted interest since the first formulation of quantum 
theory. Despite the fact that the shape of an autoionizing resonance is stipulated by quantum 



interference l20 



21| , for a long time the continuum was regarded as an incoherent dissipation 

which established that the mixing 



medium. It was until the work 0, H Q, Q, Q 
of discrete and continuous states in a strong, driving electromagnetic field led to similar 
constructive and destructive interference processes. Thus the opportunity was shown to 
produce an autoionizing-like resonance embedded in an otherwise unstructured continuum, 
which displays itself in the variety of p hoto- and nonlinear-optical processes similar to real 
autoionizing states 

HQ 0,0 HQ- 

The attractive advantage of this option is that 
such laser-induced continuum structures (LICS) can be produced, at least in principal, in any 
necessary area of a continuum, and their strength is controlled with the dressing laser. Since 



the first successful experiments 
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32l |33|, which confirmed the principal theoretical 



predictions, a great progress has been achieved in develo ping a p rofound understandin g of 



continuum coherence and relatec 
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laser-induced processes jS, 
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The appearance and consequences of coherence and interference processes are different in 
the continuous wave regime, where the relaxati on p rocesses play a crucial role, and in the 



pulsed regime of rapid adiabatic passage 
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j. Extensive studies of population transfer between two discrete states 
via the upper- lying continuum in Raman-like A-schemes have revealed that, despite the 
detrimental effects, almost complete population transfer and corresponding photoionization 
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suppression 
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is possible in such schemes, provided that the driving laser pulses are properly 
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871 l8q . It was also found that coherence processes play 



an important role in two-photon dissociation mediated by an intermediate resonance with a 



discrete level allowing for the control of photodissociation and photoassociation 
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90, 
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94| . Population transfer and dissociation play an important role as the accompanying 



processes in the three-pulse four-wave mixing techniqu e, which is us ed for spectroscopy and 



coherent quantum control in chemistry 



101|. 



96, gJQflllQO, 

Based on the outlined achievements, this paper further develops the theory of coherent 
quantum control, which employs LICS and rapid adiabatic passage. We propose and inves- 
tigate a more complex three-pulse scheme (Fig. which allows additional flexibility and 
means for such control. Unlike the so far investigated schemes, we report on novel opportu- 
nities for manipulating dissociation and population redistributions between the upper bound 
states n and / through the interplay of two LICS attributed to lower-lying energy continuum. 
Only one discrete state here is connected to the ground state by one-photon transitions. All 
radiations are assumed to be Gaussian pulses with variable sequence and duration, strong 
enough to drive molecular transitions. Corresponding equations for slowly-varying proba- 
bility amplitude are derived, and extensive numerical simulations are performed aimed at a 
demonstration of the outlined opportunities under typical experimental conditions. 

This paper is organized as follows. The set of equations for slowly-varying probability 
amplitudes pertinent to the problem under consideration are derived in Sec. ini By tur ning 
off the first or third laser, our scheme can be easily reduce to those investigated in 0, IsJ 
where excellent agreement between theory and experiment was reported. Corresponding 
analytical solutions are found for further analysis of time-dependent destructive and con- 
structive interference and for testing the theory by comparing with the known results. They 
are also useful for determining the dependence of such processes on the complex effective 
Fano parameters, which represent the properties of the specific continua. The results of var- 
ious numerical experiments are presented in in Sec. IIIIl Suppression of photodissociation 
through continuum coherence and overlap of two LICS is simulated in Sec. IIII Al The fea- 
tures of photodissociation controlled by two ordered pulses are simulated in Sec. IIIIBI We 
investigate the possibility of manipulating the dissociation yield and populations of excited 
states, driven by three pulses with variable sequence and duration. Our goal is to demon- 
strate the feasibility of achieving almost any necessary ratio between the branching yields 
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by the end of the pulses due to the interference of quantum pathways through a variety of 
continuum states. This is done with the aid of numerical experiments in Sec. IIII CI We 
show how the pulse parameters must be adjusted in order to enable such opportunities. 

II. BASIC EQUATIONS 

The proposed coupling scheme is illustrated in the energy level diagram depicted in Fig. ^ 
Radiation at frequency ui couples the bound-bound transition m — n, and radiations at uj2 
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FIG. 1: Energy levels and driving fields. 

and cus couple the bound states n and / with states of the dissociation continuum e, as shown 
in the picture. Initially, only the lowest level is assumed populated. At ui tuned near the one- 
photon resonance, stepwise and two-photon transitions interfere. The two-photon transition 
rate between the levels n and / depends on the detuning Qnf- Actually, all detunings are 
intensity-dependent, which determine the complex dynamics of the interference between 
various stepwise and multiphoton processes. We assume that each laser drives only one 
transition and accounts only for resonant coupling depicted in Fig. ^ Therefore, we ignore 
other off-resonant and incoherent processes caused by the driving field. We also assume a 
single continuum, i.e., no degeneracy, in order to avoid unnecessary complexity related to the 
processes, which are not the subject of this paper. Basically, the set of equations depends on 
whether the system is open or closed. In the case of a closed model, the lowest level m is the 
ground state, and the sum of the level populations is always equal to 1. Alternatively, in an 
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open scheme, level m is an excited state that relaxes, and the sum of the populations depletes 
with time. For the purposes of discussion, we will consider the second case, although for 
pulses considerably shorter than all the relaxation times, the difference between the closed 
and open schemes disappears. In the rotating-wave approximation, the equations for slowly- 
varying probability amplitudes in the case of an open energy-level system are written as 

ide = G£nexp(ifie„t)a„ + Gef exp{iQsft)af 
idf + i'jfaf = J dEGf£exTp{i^lfet)as, 

Here G^n = E2den/2h, G^j = E^dej/lh, and Gmn = are coupling Rabi frequencies; 

VLmn = uji — ujmn, ^ne = UJ2 — Une and Qfs = UJ3 — Ufs are frcqucucy resonance detunings; 
Uis = {Ei — e)/h; Ei and e are the energies of the corresponding discrete and continuum 
states; and is the decay rate of level i. Following [2^, we introduce the values 

af(t) = afexp{iQfsot), am{t) = amexp{iQrnsot), (2) 

where is the amplitude, which slowly varies over the continuum in the vicinity of the 
energy Eq = Ef - hjs, 

^feo — ^feo ~ ^^3 = 0, ^neo — ^neo — ^2 = ^feo ~ ^nf — ^2 = ^nfj 

^meo — ^meo — ^1 + = ^mn " ^neo — UJi + iU2 = ^mn ~ ^nf- (3) 

The energy is chosen in the vicinity of En — fiuj2, if E^ = 0. We assume that the 
coupling parameters (energy density of the transition oscillator strengths fie) vary around 
£0 substantially only over the continuum energy intervals e, which are much greater than 
the maximum of the characteristic widths of the discrete states, including power-broadening 
h\Vjk\ and pulse spectral width r~^: 

e » /i7„ h/r, h\V,u\d''f/del^,^ « 9//9£,=,„. (4) 

Such requirements are fulfilled practically for all realistic atomic and molecular continua in 
an energy range that is well above the ionization or dissociation threshold. Since interference 



processes appear at |l^,nn|, and \Qmn — ^nf\ less than or on the order of the same value, 

this enables us to separate the slow and fast variables. With the aid of (j21), we obtain from 



iae-^leeo(^e = GenCln + Gefaf, (5) 

«„-s„(Si„;-nJ = G„„ir„ + /d.G„3.. (6) 

iaf + iaf'jf = J deGfette, (7) 

iOim (^miy^mn ^nf i'~im) Gmn^^n- (8) 

The solution of © is 

Taking slowly- varying values Geitti (as compared with the oscillating exponents) out of the 
integral, we obtain 

It follows from the initial conditions that Ci = 0, and therefore the solution of (0) is 

as = -{Genan + G efU f) / Vteeo' (9) 



This can be derived also directly from (jSJ because the inequality \ae\ <^ |f2eeoa£| is correct 
for the major continuum interval. Further, with the aid of the (^-function, 

[i{ujk-ujje)Y^ = T^5{'^k-i^je)-i'P{uJk-i^je)~'^] k = 1,2,?,] j = m,n,f, (10) 

where V stands for the principal value of an integral, Eqs. (jHl)-® can be presented in the 
form 

dcirn/ dt iG^fiCln (^m ~l~ ^^meo)'^"^' 

dan/dt = -iG*^^am. - 7n./(l + iqnf)af - [7n + Inn + i{^neo + Snn)]an, (H) 

dttf/dt = -7/n(l + iqfn)an - {-if + 7// + ^5ff)af- 
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where 



= ftp / & . G„G„/(.„ - .) + Im (G„.G.„/p.,) ; 

7// = ^fiGfsoGeof + Re (GfkGkf/pkf) , 

^ff = ^'P j d^- GfeGef/{eo -e) + lm (GfkGkf/pkf) ; 

7„/ = nhGneoGeof + Re (GnkGkf/pkf) , (12) 

(^n/ = hV J de- GneGefl {eq - e) + Im (GnkGkf/Pkf) ; 

7/„ = TlhGfsoGeon + Re (GfkGkn/Pkf) , 

Sfn = hV j de- GfeGen/{eo - + Im (GfkGkn/Pkf) ; 

Pfc/ = Tfc/ + z(cjfc/ - u;3), = 5.y/7ij, j = n, f. (13) 

Besides the continuum states, the contribution of other non-resonant levels k is taken into 
account, and a sum over the repeating k index is assumed. Contributions from these levels 
may occur comparable to those of the continuum states. As seen from the equations 
the values 7jj and 6ij describe light-induced broadening and shifts of discrete resonances 
stipulated by the induced transitions between them through the continuum. The magnitude 
of the shift and its sign are determined by the overall counterbalance of the continuum states 
below and above Eq. The parameters qij = Sij/'-fij are analogous to the Fano parameters for 
autoionizing states. They characterize the relative integrated contribution of all off-resonant 
quantum states compared to the resonant ones. Within the validity of (j3)), their dependence 
on the field intensities can be neglected. The ratio of the real and imaginary parts (relative 
phase of the corresponding induced atomic oscillations) plays a crucial role in whether the 
interference is constructive or destructive. This depends on the effective Fano parameters qtj, 
which are determined by the distribution of the oscillator strengths over the continuum and 
discrete states and on the positions of the resonant continuum states, and can be controlled 
through multiphoton detunings. 

Further, we will use dimensionless variables, scaled to the pulse l-Eip half-duration r at 

the 1/e level: T = t/r, 'JiT = rji, g^n = GmnT, Onf = InfT, Qnn = Innr, Off = IffT, 
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= QmnT, ^nf = ^nf'T- Then the equations take the form 

(Mn/dT = -ig^^a„^ - 5'„/(l + iqnf)af - [rjn + gnn + i{^nf + qnn9nn)]an, (14) 
dttf/dr' = -gfn{l + iqfn)an - ijlf + 9ff + iqff9ff)af- 

With the aid of these equations, the dissociation probabihty is calculated as 

W = 2 j dT {(5(„„|a„p + 5(//|a/p + 2Re [5(„/a„a}exp(iA„/T)] } . (15) 

The first two terms, proportional to the squared moduli of the probability amplitudes, de- 
scribe the multistep dissociation associated with the populations of levels n and /. The 
third term describes quantum control through the interference of coherent quantum path- 
ways. As seen from Eqs. and ()14j) . the magnitude and sign of the interference term in 
ffT3j) strongly depends on the phases of the probabihty amplitudes and, consequently, on the 
parameter gj„. 



III. NUMERICAL SIMULATION OF TWO- AND THREE-PULSE COHERENT 
CONTROL 

Many experiments on the coherent control of branching chemical reactions have been 
carried out with sodium dimers Na2. In this case, level m of our model (Fig. can be 
attributed to the state X^S+(t; = 28, J = 10), and levels n and / to the states A^Il^{v = 
37, J = 11) and A^T,^{v = 47, J = 11), which are coupled by strong transition dipoles 
with the dissociation continuum Na(3s)+Na(3s) j^. Alternatively, vibrational states of the 
electronic excitation B^Ilu{v, J = 11) can be chosen as level /. The characteristic relaxation 
rates of these states are: 7^ = 2 • 10'' c~^, 'jn = If = 1-2 ■ 10® c~^. We assume all pulses to 
be Gaussian, 

9mn = ^Lexp(-TV2), 

7nn = 7°nexp[-(r-A2)V4], (16) 

7// = 7?/exp [-{T-Asr/dl], 

where A2 and A3 are the delays of the pulses E2 and with respect to the pulse Ei, 
<^2 = T2/T and d-s = t^/t are the durations of pulses E2 and E^ {t2 and Ts) scaled to the 

8 



duration of the first pulse r. These values must be selected within a time domain range 
shorter than r = 10~^s in order to avoid relaxation over the whole period of excitation. 
As seen from Eq. (jl5p . the overall analysis of the control reduces to an analysis of the 
interference and multistep terms. 



(a) (b) 




FIG. 2: Dissociation (a) and population transfer from level n to level / (b) vs time and delay 
of the pulse E2 relative to (A2 = 0, A3 is negative ). The field Ei is turned off, the initial 
population of level n equals 1, = 3.61, g^j = 9.61, g^j: = 5.89, the pulse durations are equal 
and much shorter than the relaxation rates, r2„j = 0, and the Fano parameters are Qnn = 0.2, 

'iff = -0-5' Inf = lO- 



A. Suppression of photodissociation by continuum coherence and overlap of two 
Lies 

Figures El and El demonstrate feasibilities of manipulating constructive and destructive 
interference aimed at control of the population transfer between the upper levels n and /. 
Here the first field is turned off and only level n is initially populated. Figure El simulates the 
population dynamics when the carrier frequencies of the fields are tuned to the two-photon 
resonance with Ufn- It shows a dramatic dependence of the output on the temporal shift 
between the pulse peaks. Such dynamics is stipulated by the time-dependent shift of the 
two-photon resonance, which determines constructive or, alternatively, destructive interfer- 
ence. In other words, the position, strength and overlap the two LICS vary in time, which 
bring about the dynamics and consequences of the interference. For the given parameters, 
we have almost complete suppression of dissociation [plot (a)] and about 80% population 
transfer to level / [plot (b)]. This occurs for counterintuitive order of pulses, when third 
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FIG. 3: Dependence of the population transfer (a,b) and dissociation (c,d) on the Fano parameter 
Qnf (a,c) and on two-photon detuning A„j (b,d). A2 = 0, A3 = —3.9 for all plots. (a,c) A„j = 0, 
(b,d) Qnf = 10. All other parameters are the same as in Fig. [3 

field reaches its maximum at about four pulse half-duration before the maximum of second 
pulse. Alternatively, almost complete photodissociation can be ensured with different pulse 
shifts. Thus our simulations demonstrate a behavior, which is in a erood agreement with the 
cclusio,. a.d ob.e„at,„,. of o.he. stud.es QQ. Fig.eE|p.ot. (a) a.d (c), d,spla,s a 
strong dependence of dissociation and population transfer on the effective Fano parameter 
Qnf- This because they determine the magnitude and sign of the resonance shift (position 
and shape of LICS). The plots (b) and (d) show that similar manipulations can be performed 
by an appropriate adjustment of the two-photon detuning. 

B. Photodissociation controlled by two ordered pulses 

Figure |3] shows the features of dissociation from the lowest state controlled with two 
pulses. Here the third field is turned off, only level m is initially populated, the carrier 
frequency of the first field is tuned to resonance with the transition mn, and dissociation is 
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FIG. 4: Dependence of the two-photon dissociation (£^3 = 0) (a,b,c) and population of the ground 
level (d) on the intensity and pulse delay of the driving fields. ^Imn = 0, = 2. (a) g^^ = 3.61; 
(b) g^^ = 400; (c,d) A2 = 0. All other parameters are the same as in Fig. [2 

caused by the interference of two-step and two-photon processes. The plots indicate that, 
due to the time-dependent shift of one-photon resonance induced by the second field, the 
dynamics of photodissociation and its dependence on the pulse sequences and shift changes 
considerably with the increase of the peak intensity of the second field [plots (a) and (b)]. 
At a given pulse-shift there exists an optimal intensity of this field [plot (c)]. The increase 
of the intensity of the second field leads to decoupling of the first one with the transition 
mn [plot (d)] and, consequently, to a decrease in the overall dissociation [plot (c)]. Rabi 
oscillations of the populations at the transition mn in the range of relatively small intensities 
of the second field [plot (d)] may give rise to substantial modification of the dynamics and 
to suppression of dissociation [plot (b)]. 
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FIG. 5: Dissociation and population of levels driven by three laser pulses, ftmn = ^nf = 0, 
= 2, g^^ = 0.25, g^j = 0.36, g^j = 0.3, d2 = 1, = 1.6; A3 = 0. The Fano parameters are 
Qnn = 0.2, Qff = —0.5, Qnf = 10. Initially, only level m is populated, (a) and (b) - dependence 
on the pulse delay, (c) and (d) - dynamics of the populations and dissociation (S|ajp is the sum 
of populations including the energy integrated population of the continuum W). (c) A2 = 0, (d) 
A2 = -1.5. 
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(a) (b) 
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FIG. 6: Dependence of the three-pulse driven population of level / on the Fano parameter g„j (a) 
and on the two-photon detuning A.„j (b). A2 = 2.8, A3 = 0. (a) 0^/ = 0, (b) Qnj = 10. All other 
parameters are the same as in the previous figure. 

(a) (b) 




FIG. 7: Dependence of the dissociation and of the population of the intermediate level n on the 
two-photon detuning, where Qnf = 10. All other parameters are the same as in the previous figure. 

C. Three-pulse control of photodissociation and population transfer by continuum 
coherence 

With an understanding of the features described above, we shall investigate the dynamics 
of more complex scheme which embraces both schemes considered above by assuming all 
three fields at c^i, 002 and (Fig. turned on. First, we shall simulate the case when one-. 
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two- and three-photon resonances are fulfilled {Qmn = ^nf = 0), and the first and second 
pulses are peaked at the same instant as the third pulse, which is somewhat longer. The 
second pulse is of the same duration as the first one and can be either advanced or delayed 
with respect to the other two. 

Figure El [plots (a) and (b)] shows that the dissociation output and population of the 
upper level can be controlled within the wide intervals by adjustment of the sequence and 
the time shift of the second pulse. With two examples [plots (c) and (d)], the simulation 
demonstrates the feasibility of achieving either comparable populations of all three discrete 
and integrated continuum states by the end of the pulses [plots (c), where all pulses peak at 
the same instant] or, alternatively, of ensuring a dominant population of the upper discrete 
state and dissociation output by advancing the second pulse regarding the other two [plot 
(d)]. Here, S|ajp is the sum of populations including the energy integrated population of the 
continuum W. Figure IHl [plot (a)] shows that the dynamics of the populations is considerably 
different for different Fano parameters (keeping the other parameters the same). It also shows 
that, for fixed time delays, the output can be controlled by the adjustment of appropriate 
two-photon detuning [plot (b)]. Figure [3 shows similar feasibilities of manipulation of the 
populations of the intermediate level and dissociation. Dark resonances, characterized by 
suppressed dissociation and selective population redistribution between the coupled bound 
states, or, alternatively, enhanced dissociation are shown to be achieved by variation of 
two-photon (and consequently three-photon) detuning. 



IV. CONCLUSION 



The theory of adiabatic passage between three discrete levels and a low-lying dissociation 
continuum is developed. Possible quantum control of two-photon dissociation (A-scheme) us- 
ing auxiliary short laser pulses at adjacent bound-free transitions (^-configuration) is shown. 
The proposed three-pulse scheme enables one to manipulate dissociation through the dark 
states not connected with the ground one by the allowed transition. Besides dissociation, 
the scheme under investigation enables one to control the population transfer between two 
upper discrete levels via the lower-energy dissociation continuum, while a direct transition 
between these states is not allowed. A system of coupled equations for the probability am- 
plitudes is employed, where three time-shifted pulses drive strongly the quantum system, so 
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that nonlinear quantum interference occurs through the overlap of two laser-induced contin- 
uum structures. The opportunities of manipulating photoinduced processes by judiciously 
frequency-detuned and time-delayed pulses are explored through extensive numerical simu- 
lations with the aid of a model relevant to typical experiments. The simulations display good 
agreement with the known experiments on two-pulse population transfer between metastable 
states via the ionization continuum and with experiments on two-photon photodissociation. 
The dependence of the dissociation spectrum on the magnitudes of the composite Fano pa- 
rameters for excited levels, on detuning between two laser-induced continuum structures, on 
intensities, and on the counterintuitive order and delay of the pulses are investigated. We 
have demonstrated the opportunities of considerable suppression of photodissociation and 
considerable transfer of the lower-state population to the excited states or, alternatively, al- 
most complete photodissociation, or redistribution the population between continuous and 
discrete states with the required ratio by the end of the pulses. The gained efficiency and 
flexibility may compensate for the increase in experimental complexity. 
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